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ABSTRACT 

We have mapped the warm molecular gas traced by the H2 S(0) — H2 S(5) pure 
rotational mid-infrared emission lines over a radial strip across the nucleus and 
disk of M51 (NGC 5194) using the Infrared Spectrograph (IRS) on the Spitzer 
Space Telescope. The six H2 lines have markedly different emission distributions. 
We obtained the H2 temperature and surface density distributions by assuming 
a two temperature model: a warm (T = 100 — 300 K) phase traced by the 
low J (S(0) - S(2)) lines and a hot phase (T = 400 - 1000 K) traced by the 
high J (S(2) - S(5)) hues. The lowest molecular gas temperatures are found 
within the spiral arms (T ~ 155 K), while the highest temperatures are found 
in the inter-arm regions (T > 700 K). The warm gas surface density reaches a 
maximum of 11 Mq pc^^ in the north-western spiral arm, whereas the hot gas 
surface density peaks at 0.24 Mq pc^^ at the nucleus. The spatial offset between 
the peaks in the warm and hot phases and the differences in the distributions 
of the H2 line emission suggest that the warm phase is mostly produced by 
UV photons in star forming regions while the hot phase is mostly produced by 
shocks or X-rays associated with nuclear activity. The warm H2 is found in the 
dust lanes of M51, spatially offset from the brightest Ha regions. The warm 
H2 is generally spatially coincident with the cold molecular gas traced by CO 
(J = 1 - 0) emission, consistent with excitation of the warm phase in dense 
photodissociation regions (PDRs). In contrast, the hot H2 is most prominent 
in the nuclear region. Here, over a 0.5 kpc radius around the nucleus of M51, 
the hot H2 coincides with [O IV] (25.89 fim) and X-ray emission indicating that 
shocks and/or X-rays are responsible for exciting this phase. 
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1. Introduction 

Star formation and galactic evolution are connected via the molecular gas in a galaxy. 
In the Milky Way, star formation occurs in molecular clouds, although not all clouds are 
actively forming stars. On a global, galactic scale, star formation may be triggered when- 
ever the molecula r gas surface density is enhanced, for example, by a spiral density wave 



(IVogel et al. 19881). by increased pressure or gas density in galactic nuclei (I Young and Scoville 1991 



Sakamoto et al. 19991: ISheth et al. 20051). bv hydrodynamic shocks along the leading edge of 



bars (ISheth et al. 2000l : ISheth et al. 2002[) , and in the transition region at the ends of bars 
( iKenny and Lord 199ll : ISheth et al. 20021 ). How does this star formation affect the surround- 
ing molecular gas? How is it heated and what is the distribution of the gas temperatures? 
How does the mass of the warm and hot gas vary from region to region? We address these 
questions using spectral line maps from a radial strip across the grand-design spiral galaxy, 
M51. 

M51 (the Whirlpool galaxy, NGC 5194) is a near by, face-on s piral galaxy that is rich in 



molecular gas. Its proximity (assumed to be 8.2 Mpc (ITuUy 19881 ) ). face-on orientation, and 
grand-design spiral morphology make it the ideal target for studies of the interstellar medium 
(ISM) across distinct dynamical, chemical, and physical environments in a galaxy. Studies of 
the molecula r gas within M51 h ave revealed giant molecu l ar associations (G MAs) along the 



spiral arms (jVogel et al. 19881 : iRand and Kulkarni 1990l : lAalto et al. 19991 ). a reservoir of 



molecul ar gas in the nuclear region that is massive enough to fuel the active galactic nucleus 



(AGN) (IScoville et al. 19981) . and spiral density wave triggered star- formation in molecular 



clouds (IVogel et al. 19881 ). In addition to being well-studied at millimeter and radio wave- 



limeter wavelengths ( 


Pa 


lumbo et al. 19851: 


, „ ■ , „^,v,^ ^, ^ ^^^^ 

Terashima et al. 1998 


Matsushita et al. 2004: 


Meiierink et al. 2005h. 



In this paper we present maps of the H2 S(0) — H2 S(5) pure rotational mid-infrared 
lines over a strip across M51 created from Spitzer Space Telescope Infrared Spectrograph 
(IRS) spectral mapping mode observations. The mid-infrared H2 lines trace the warm 



temperature, mass 


(RieoDoulou et al. 2002: 


Higdon et al. 2006 


(Neufeld et al. 1998: 


lNe„feld et al. 2nnfi) across th. M51 strins 



and ortho-to-para ratio 
We use the inferred dis- 
tributions to place constraints on the energy injection mechanisms (i.e. radiative heating, 
shocks, turbulence) that heat the warm molecular gas phase of the ISM. 



"'^ While we are always exploring the warm phase of H2, in this paper we refer to a warm and a hot phase 
corresponding to temperatures of T = 100 - 300 K and T = 400 - 1000 K, respectively. 
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2. Observations and Data Reduction 
2.1. Spectral Data 

We mapped a radial strip across M51 using the short-low (SL; 5 - 14 .5 /xm) and long-low 



(LL; 14 - 38 /im) modules of the Spitzer IRS in spectral mapping mode (iHouck et al. 20041 ) . 
The radial strips were 324" x 57" and 295" x 51" in the SL and LL, respectively. Each 
slit position was mapped twice with half-slit spacings. In total, 1,412 spectra were taken in 
the SL and 100 were taken in the LL. Integration times for individual spectra were 14.6 s 
in both the SL and LL. Dedicated off-source background observations were taken for the SL 
observations. Backgrounds for the LL observations were taken from outrigger data collected 
while the spacecraft was mapping in the adjacent module. Figured] presents the astronomical 
observation requests (AORs) overlaid on the Spitzer Infrared Array Camera (IRAC) 8 /xm 
image of M51. 

The spectra were assembled fror n the basic calibration data (BCD) into spectral data 



cubes for each module using CUBISM ( iKennicutt et al. 20031 : ISmith et al. 2004 



Smith et al. 2007af ) 



Background subtraction and bad pixel removal were done within CUBISM. The individual 
BCDs were processed using the S14.0 version of the Spitzer Science Center (SSC) pipeline. 
In CUBISM, the SL and LL data cubes have 1'.'85 and 5"08 pixels, respectively. The pixel 
size is half the full width at half maximum (FWHM) of the point spread function (PSF) 
at the red end of a given module. In principle, the PSF should vary with wavelength but 
since the PSF is undersampled at the blue end of the module, it is approximately constant 
across a given module. So the approximate resolution of the SL and LL modules and maps 
of spectral features observed in the SL and LL modules is 3'.'7 and 10"1, respectively. 

We created continuu m-subtracted line fl ux maps of the H2 S(0) — H2 S(5) lines using a 



combination of PAHFIT (ISmith et al. 2007bl ) and our own code. PAHFIT is a spectral fltting 
routine that decomposes IRS low resolution spectra into broad PAH features, unresolved line 
emission, and grain continuum with the main advantage being that it allows one to recover 
the full line flux of any blended features. Several H2 lines are blended with PAH features 
and atomic lines in IRS low resolution spectra: H2 S(l) with the 17.0 /im PAH complex, 
H2 S(2) with the 12.0 and 12.6 /im PAH complexes, H2 S(4) with the 7.8 and 8.6 /im PAH 
complexes, and H2 S(5) with the [Ar II] (6.9 /im) line. PAHFIT also solves for the foreground 
dust emission and dereddens the emitted line intensities. Our code concatenates SLl and 
SL2, and LLl and LL2 data cubes into two cubes, one for SL and one for LL and smoothes 
each map in the cubes by a 3 x 3 pixel box, conserving the flux, to increase the signal-to- 
noise ratio of the spectra. Then, for each pixel, the spectrum is extracted and PAHFIT is 
run to decompose it. Our code saves the location of the pixel on the sky along with the 
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PAHFIT output (i.e. integrated line flux, line FWHM, line equivalent width, the uncertainty 
in the line flux, the flt to the entire spectrum and the flt to the continuum) for each spectrum 
and uses this information to construct line flux maps for all of the mid-infrared features. In 
addition to creating line maps, our code creates maps of line FWHM, line equivalent width, 
and uncertainty in the flux and data cubes of the flt to the entire spectrum, the flt to the 
continuum, a continuum-subtracted data cube, and a residual data cube. These parameters 
(the uncertainty in the line flux, line FWHM, and line equivalent width) are all determined 
within PAHFIT. For a more thorough discussion of how PAHFIT determines quantities such 
as the uncertainty in the integrated line flux, line FWHM, and line equivalent width, see 
Smith et al. (2007b). 

Figure [2] presents sample Spitzer IRS low resolution spectra and the PAHFIT spectral 
decomposition for 3 different single pixel regions across the M51 strip. The SL and LL 
spectra are plotted separately because we decomposed the SL and LL data cubes separately. 
The locations of the regions are marked on the maps of the H2 surface density in Figure [51 
While we have decomposed the entire spectrum and created maps of all of PAH features and 
spectral lines, in this paper we focus primarily on the H2 line maps. 



2.2. Ancillary Data: CO (J = 1 - 0), Optical V Band, Ha, and X-ray 

Observations 

In this section we briefly discuss the ancillary data that we have used in order to un- 
derstand H2 excitation in M51. We take note of the image resolutions for comparison to 
the Spitzer IRS beam. The Berkely Illinois Maryland Array (BIMA) C O (J = 1 - 0) map 



was acquired as p art of the BIMA Survey of Nearby Galaxies (SONG) ( iRegan et al. 2001 



Heifer et al. 2003f ). At the distance of M51, the SONG beam (5'.'8 x 5'.'1) subtends 220 pc x 
190 pc. The optical V band and continuum-subtracted Ha -|- [N II] images of M51 were ob- 
tained from the Spitzer Infrared Nearby Galaxies Survey (SINGS) archiveEl The native 
pixel scale for both images is O'.'S and the angular resolutions are ~ 1". X-ray emission from 
M51 was observed by the Advanced CCD Imaging Spectrometer (ACIS) on the Chandra 
X-Ray Observatory on 20 June 2000. The resolution of the image is ~ 1". The X-ray image 
that we use has been presented and discussed in Terashima and Wilson (2001). 



^http: / /data. spitzer. caltech.edu/popular/sings/ 
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3. Results 

3.1. The Distribution of H2 Emission 

We have detected and mapped H2 emission from the six lowest pure rotational H2 lines 
(Figure [3]). The maps reveal remarkable differences in the distribution of H2 emission in 
M51. H2 S(0) emission is strongest in the northwest inner spiral arm peaking at a flux of 
3.7 X 10^^^ W m~^ and decreases by a factor of 2 in the nuclear region. In contrast, the H2 
S(l) emission peaks in the nucleus of the galaxy at 1.0 x 10~^^ W m~^ and has an extension 
of equal brightness towards the northwest inner spiral arm. In the spiral arm itself, the H2 
S(0) peak is offset from the H2 S(l) emission by 10" (~ 380 pc). We detect emission from 
the H2 S(0) and H2 S(l) lines to the outer limit of our radial strip, ~ 6 kpc from the nucleus 
of the galaxy. In the outer spiral arms, the H2 S(0) flux is a factor of 2 times lower than in 
the inner northwest spiral arm and the H2 S(l) flux is a factor of 5 times lower than in the 
nucleus. 

The H2 S(2) — H2 S(5) maps show different molecular gas distributions within M51 
through each H2 line. The brightest H2 S(2) emission is from the nucleus at 2.2 x 10"^^ W 
m~^. We also see bright H2 S(2) emission from the inner northwest spiral arm at half the 
flux of the nuclear peak. The H2 S(3) peak at the nucleus is 1.4 x 10^^"^ W m~^, a factor 
of ~ 6 greater than the H2 S(2) nuclear peak. There is also a linear bar-like structure in H2 
S(3) emission across the nucleus of the galaxy at a PA^-IO*^ and length of ~ 1.5 kpc. The 
emission peaks in the H2 S(2) and H2 S(3) maps are not spatially coincident. For instance, in 
the northwest inner spiral arm the brightest H2 S(3) emission is further down the spiral arm 
than the brightest H2 S(2) emission. Offsets like these suggest variations in the excitation 
temperature from region to region within a galaxy, and even within a spiral arm. 

The H2 S(4) and H2 S(5) lines are brightest at the nucleus with fluxes of 3.1 and 8.0 
X 10~^^ W m~^, respectively. The H2 S(4) line shows emission in the nucleus and in the 
spiral arm to the west. In the spiral arm to the west of the nucleus, the H2 S(4) flux is 2.1 
X 10^^^ W m~^. This is notable because the spiral arm to the southwest of the nucleus is 
very bright in CO and and studi es have revealed very high molecular gas column densities in 



the southwest inner spiral arm (ILord and Young 1990l : lAalto et al. 19991 ). H2 S(5) emission 



is asymmetric in the nucleus and mimics the morphology of the H2 S(3) line with extended 
emission to the north of the nucleus. The trail of intensity peaks along the long northeast 
edge of the H2 S(4) and H2 S(5) maps is not real. It is likely due to difficulty smoothing the 
maps near the edges. These regions do not effect the determination of the H2 temperature 
and surface density due to the offset of the SL strip relative to the LL strip. 



Based on the uncertainty maps generated while creating the H2 line maps, errors in the 
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H2 S(0) line flux that range from 15% to 75% with the largest uncertainties being found in 
the inter-arm regions and the outer northwest and southeast spiral arms. Errors in the H2 
S(l) line fluxes range from 10% to 70% with the largest uncertainties being found at the 
distant northwest and southeast portions of the M51 strip. Errors in the H2 S(2) and H2 
S(3) line fluxes range from 10% to 85% with the largest uncertainties being found toward 
the inter-arm regions. Errors in the H2 S(4) and H2 S(5) line fluxes are ~ 30% and 20% 
within a 1 kpc radius of the nucleus and become higher moving radially away. 



3.2. Mapping H2 Excitation Temperature and Surface Density across M51 

3.2.1. Modeling H2 Excitation Temperature and Surface Density 

The pure rotational lines of molecular hydrogen provide a powerful probe of the condi- 
tions of the ISM by placing constraints on the energy injection that excites H2. For example, 
Neufeld et al. (2006) discuss shock excitation of H2 and Kaufman et al. (2006) discuss H2 
excitation in photodissociation regions (PDRs). Using the maps of H2 emission, we modeled 
the H2 temperature and mass distribution over the radial strip in M51 following the methods 
described in Rigopoulou et al. (2002) and Higdon et al. (2006). 

First, we smoothed the H2 S(l) — H2 S(5) maps to the resolution of the H2 S(0) map, 
lO'.'l. The maps were then interpolated to the same spatial grid. Excitation diagrams across 
the strip were derived from the Boltzmann equation 

N,/N = ((7(z)/Z(Tex))exp(-T,/Tex) (1) 

where g{i) is the statistical weight of state i, Z(Tex) is the partition function, Ti is the energy 
level of a given state, and Tex is the excitation temperature. N and Nj are the total column 
density and the column density of a given state i and Ni is determined directly from the 
measured extinction-corrected flux by 

iVi = 47r X flux{{}/{nA{{}hiy{i)) (2) 

where A{i) is the Einstein A-coefficient, z/(i) is the frequency of state i, Q is the solid angle of 
the beam, and h is Planck's constant. Table [1] lists the values for the wavelength, rotational 
state, Einstein A-coefficient, energy, and statistical weight of the pure rotational levels of 
H2. 

In order to derive temperature and surface density distributions we assume a two tem- 
perature model for the H2. To determine the hot (T = 400 - 1000 K) phase temperature, 
we do a least squares fit to the H2 S(2) - II2 S(5) column densities in the excitation diagram 
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at every pixel in our maps. We then subtract the contribution of the hot phase from the 
lower J lines and do a least squares fit to the column densities of H2 S(0) - H2 S(2) lines at 
every pixel in our maps to determine the temperature distribution of the warm (T = 100 - 
300 K) phase. The warm and hot phase surface density distributions are derived from the 
column densities of the H2 S(0) and H2 S(3) lines, respectively. The column densities are 
easily converted to mass surface density by determining the mass of the warm and hot H2 
within every pixel. 

Figure H] shows excitation diagrams and the fits to the warm and hot H2 phases for 
three different regions across the M51 strip. The three regions are marked on the H2 sur- 
face density maps in Figure [51 In the nuclear region (Region 2) the ortho and para levels 
appear to lie along the same curve indicating an ortho-to-para ratio (OPR) of 3; the ex- 
citation diagrams do not exhibit the "z igzag" characteristic of a non-equlibrium H2 OPR 
( jNeufeld et al. 19981 : iFuente et al. 19991 ). The slope of the curve appears to decrease as the 
rotational state increases indicating that we are sampling a continuous range of H2 temper- 
atures with the higher rotational states (S(2) - S(5)) probing a hotter phase of H2 than the 
lower rotational states (S(0) - S(2)). In the southeast and northwest spiral arms (Regions 1 
and 3, respectively) the lower J (II2 S(0) - H2 S(3)) levels exhibit an OPR of ~ 3. The H2 
S(4) measurement shows significant scatter in the excitation diagrams outside of the nuclear 
region of M51. This would indicate that the OPR is less than 3, however, due to the low 
signal-to-noise ratio of the H2 S(4) map, we do not believe that the OPR determined from 
the H2 S(4) fiux refiects the OPR of the warm II2. We assume an OPR of 3 when deriving 
the H2 temperature and mass distributions, which is consistent with Roussel et al. (2007) 
who found an OPR of 3 for the nuclear region and four of seven HII regions studied in M51. 



3.2.2. Warm and Hot H2 Excitation Temperature and Surface Density Distributions 

Figure [5] presents the warm [left) and hot {right) H2 surface density distributions 
across the M51 stripl^l The highest gas surface density for the warm H2 phase is in the inner 
northwest spiral arm at 11 M0 pc~^. The gas surface density in the outer northwest and 
southeast spiral arms is maximum at the center of the spiral arms at 3.5 Mq pc-2 and 1.0 
Mq pc~^ respectively. The hot phase surface density is highest in the nucleus and interior 
to the inner spiral arm at 0.24 Mq pc~^. The gas surface density of the hot phase in the 
spiral arms is 3 - 5 times lower than that of the nuclear region. 

Figures [6] and [7] show the distribution of the warm and hot H2 temperatures (in grayscale) 



■^Note that the non-rectangular shape of the strip is due to the offset between the SL and LL strips. 
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with the contours of the warm and hot H2 surface densities overlaid, respectively. In both 
cases we see that the temperature and surface density are inversely correlated with the 
hottest temperatures corresponding to regions of lowest surface density. For both the warm 
and hot phases, we see the temperature is higher in the inter-arm regions than in the spiral 
arms. This is real and does not result from lower signal-to-noise in the inter-arm regions. 
We tested this by extracting spectra over 0.76 kpc^ inter-arm regions between the northwest 
inner arm and northwest arm and between the southeast inner arm and southeast arm and 
found warm H2 temperatures of 177 and 175 K and hot H2 temperatures of 691 and 690K, 
respectively. These temperatures are higher than the 150 - 165 K warm H2 temperatures 
and the 550 - 650 K hot H2 temperatures found in the spiral arms. 

Figure [S] compares the warm (in grayscale) and hot (in contours) H2 surface density 
distributions. The warm H2 mass distribution peaks in the northwest inner spiral arm and 
the hot H2 mass distribution peaks in the nucleus and in the region interior to the northwest 
spiral arm. The warm-to-hot H2 mass ratio is not constant across the galaxy but is lowest (~ 
12) in the nucleus of the galaxy and increases to 170 and 136 in the southeast and northwest 
spiral arms, respectively. 



4. Discussion 

4.1. Effects of Beam Averaging on Studies of Extragalactic H2 

Previous studies have used apertu re-averages over entire galactic nuclei to derive the 



physical conditions of the molecular gas (IRigopoulou et al. 2002 : iHigdon et al. 20061 : iRoussel et al. 20071 ) 



In M51, Roussel et al. (2007) find that within the central 330 arcsec^ (4.61 x 10^ pc^), the 
warm H2 phase has a temperature of 180 K and a surface density of 3.2 Mq pc^^ (a total 
mass of M^arm = 1-5 X 10^ M0). Roussel et al. (2007) have also measured the temperature 
of the hot phase (though they do not measure the mass in the hot phase) and find a hot H2 
temperature of 521 K. 

Having spatially resolved spectra over a strip across M51, we can investigate the behavior 
of the warm and hot H2 phases on smaller scales. In the nuclear region of M51 we see that 
the warm phase temperature peaks at 192 K and decreases radially towards the inner spiral 
arms. The warm H2 surface density at the nucleus is 4.4 Mq pc~^ and decreases over a 0.5 
kpc radius surrounding the nucleus. To check the consistency of our results against those of 
Roussel et al. (2007), we averaged the warm phase temperature over a similar 412 arcsec^ 
(5.76 X 10^ pc^) aperture and found that the warm phase temperature and surface density 
are 186 K and 2.8 pc~^, respectively. Our results are consistent with previous studies 
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of warm H2 done at lower resolution and suggest that previous studies of the warm H2 
temperature and mass have yielded average values rather than a maximum (or minimum) 
values. 



4.2. Distinguishing the H2 Excitation Mechanisms 

The warm-to- hot H2 surface density ratio varies across M51 suggesting that the H2 ex- 
citation mechanisms have different effects on the warm and hot phases. The largest warm 
H2 surface densities are found in the spiral arms suggesting that the warm phase is asso- 
ciated with star formation activity. The largest hot H2 surface densities are found in the 
nuclear region of M51. This suggests that where nuclear activity is the dominant excitation 
mechanism, a larger fraction of the H2 can be maintained in the hot phase and that H2 is 
generally more efficiently heated by nuclear activity. This is further supported by the fact 
that the hot H2 surface density is ~ 660 K in the nuclear region and decreases to ~ 530 
K and 550 K in the inner northwest and southeast spiral arms, respectively. Roussel et al. 
(2007) find that H2 is generally heated by massive stars in PDRs, however, Seyferts and 
LINERs show evidence for the dominance of other excitation mechanisms such as X-rays or 
shocks. By comparing the spatial distribution of the warm and hot H2 to optical V band 
imagery. Ha, and CO (J = 1 - 0) emission, we can investigate the relationship between star 
formation and H2. Additionally, by comparing the warm and hot H2 distributions to X-ray 
observations and the spatial distribution of [O IV] (25.89 /im) emission, we can investigate 
X-ray and shock heating of H2 within M51. 

Kaufman et al. (2006) show that within galaxies, where the telescope beam size is 
generally kiloparsecs across, H2 emission could serve to probe the average physical conditions 
in the surfaces of molecular clouds. In Figure [9], we compare the warm [left) and hot II2 
[right) mass distributions to the cold molecular gas traced by CO (J = 1 — 0) emission. 
The warm and hot H2 phases appear to trace the bright CO emission in the northwest and 
southeast spiral arms. Comparison of CO to the individual H2 S(0) — II2 S(3) line intensity 
maps in Figure [10] also shows that the H2 in the spiral arms traces the bright CO emission. 
In the spiral arms, both the warm and hot H2 phases are found in PDRs and are associated 
with star formation activity. The most striking result is that in the inner spiral arms, we see 
that the CO is offset toward the nucleus from the warm H2 mass (in Figure [H]). The offset 
between the peaks in CO and warm H2 mass is 7" in the northwest inner spiral arm and 5" 
in the southeast inner spiral arm. We believe that these offsets are real with one possible 
explanation being that the H2 is tracing the regions of active star-formation within the giant 
molecular associations. 
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In Figures [TT] and [T2l we compare the warm [left) and hot [right) H2 mass distributions 
to Ha emission and an optical V band image, respectively . The Ha image reveals the 
brightest HII regions and the V band image shows the dust lanes in M51. In general, the 
warm and hot H2 concentrations are not cospatial with the brightest Ha emission regions in 
the spiral arms with the one exception being that the warm H2 mass in the northwest and 
southeast inner spiral arms appears to trace the Ha emission. The warm H2 mass contours 
show that local peaks in H2 mass are found within the dust lanes, with the Ha emission 
spatially offset towards the leading edge of the spiral arms. An example of this is in the 
northwest spiral arms where we see the H2 mass offset from the Ha spiral arms with local 
peaks being found in the dust lanes. In Figures [13] and [HI we compare the H2 S(0) — H2 S(3) 
line intensity maps to Ha emission and V band imagery, respectively. Comparison of the 
H2 S(0) map to Ha reveals that the strongest H2 emission in the northwest and southeast 
inner spiral arms is coincident with Ha emission; however, the other H2 S(0) spiral arms 
show the strongest emission in the dust lanes, offset from the Ha spiral arms. The largest 
offsets are seen in the southeast spiral arm where the H2 S(0) emission is offset from the Ha 
spiral arm by ~ 15" (560 pc). H2 S(l) emission appears to follow the dust lanes and the H2 
S(l) intensity subsides into the Ha spiral arms. H2 S(2) and H2 S(3) emission is also found 
in the dust lanes; however, there are instances (such as in the southeast spiral arm) where 
the H2 emission appears to be found straddling the dust lane and Ha spiral arm. 



The [O IV] (25.89 fim) line can be e xcited in fast shocks (v. > 100 km s"^) fiLutz et al. 19981) 



by photoionization in Wolf-Rayet stars (ISchaerer and Stasinska 19991 ). or by an AGN JSmith et al. 2004h . 



We have mapped the blended [O IV] (25.89 /xm) + [Fe II] (25.99 /im) line across M51. While 
the map contains the contribution of both lines, we have examined the SINGS high resolu- 
tion spectra and find that averaged over a 365 arcsec^ region around the nucleus of M51, 
the [O IV] line flux is 2.7 times higher than the [Fe II] fiux0 In the high resolution data 
cubes of the nuclear region, the [O IV] line flux dominates the [O IV] + [Fe II] blend and 
the [O IV] emission distribution derived from the low resolution spectra accurately reflects 
the distribution seen in the high resolution spectra. In Figure [151 we compare the [O IV] 
surface brightness to the warm (left) and hot (right) H2 distributions. The [O IV] emission 
is brightest in the nuclear region at 1.2 x 10^^^ W m~^ and within a 0.5 kpc radius of the 
nucleus, the peak is coincident with the peak in the mass of the hot H2. [O IV] surface 
brightness decreases from the nucleus to the inner spiral arm by 50 %. We resolve weaker 
[O IV] emission within the warm and hot H2 spiral arms. The [O IV] surface brightness is 
a factor of ~ 6 lower in the spiral arms than the peak intensity found in the nucleus. The 
[O IV] emission within the nuclear region of M51 is likely due to the weak Seyfert 2 nucleus 



SINGS data cubes and spectra for M51 can be found at http://irsa.ipac.caltech.edu/data/SPITZER/SINGS/galaxies/ngcl 
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dPord et al. IQSsh and is possibly associated with shocked gas from the outflows of the AGN. 



The peak of the [O IV] emission coincides with the nuclear peak in hot H2 mass, indicating 
that the hot H2 phase in the nuclear region of the galaxy is AGN or shock heated. In the 
nuclear region we observe a factor of 12 times greater warm H2 mass than the hot H2 mass. 
The warm H2 mass is much greater within the spiral arms than within the nucleus and the 
warm-to- hot mass ratio is lowest in the nuclear region where the [O IV] flux is greatest. In 
general, nuclear activity and shocks heat the ISM more efficiently causing the higher surface 
density of hot H2. 

X-ray studies of M51 have revealed bright X-ray emission from t he nucleus, the extranu- 



clear cloud (XNC, to the south of the nucleus) , and the northern loop (ITerashima and Wilson 2001 



Maddox et al. 20071 ). A radio jet that is believed to be shock heating the ISM has been 



observed emanating from the south of the nucleus towards the XNC in 6 cm imagery 



( 1 Crane and van der Hulst 19921 ). In Figure [T6| we compare the smoothed 0.5 - 10 keV band 
X-ray image to the warm (left) and hot (right) B.2 mass distributions. The 0.5 - 10 keV 
band has been smoothed to the resolution of the warm and hot II2 mass distributions and 
the nucleus, XNC, and northern loop are indistinguishable in the smoothed image. X-ray 
emission is brightest in the nucleus and decreases into the northwest spiral arm that contains 
the greatest H2 mass. There appears to be very little connection between the 0.5 - 10 keV 
X-ray band and the warm H2 mass distribution. The peak in X-ray emission is coincident 
with the hot H2 mass peak. The brightest 0.5 - 10 keV X-ray emission originates from the 
nucleus and is oriented north-to-south, similar to the [O IV] (25.89 fim) emission. The peak 
in X-ray emission is located within the peak in hot H2 mass suggesting that X-rays play an 
important role in producing the hot H2 phase. While there is a correlation between X-ray 
emission and the hot H2 phase, H2 excitation by X-rays cannot be distinguished from H2 
excitation by shocks. In Figure [T71 we compare the X-ray surface brightness to the H2 S(2) 
— H2 S(5) maps and find that the nuclear H2 emission appears to be correlated with the 
X-ray source. 



5. Conclusions 

We have spectrally mapped a strip across M51 using the Spitzer IRS low resolution 
modules. We used the spatially resolved spectra to map H2 S(0) — H2 S(5) lines across the 
strip. We find: 

1. The distribution of H2 emission in M51 varies with II2 rotational level. II2 S(0) emission 
is brightest in the spiral arms of the galaxy while the higher J transitions show the strongest 
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emission towards the nucleus. The H2 S(l) hue is brightest in the nuclear region and is offset 
from the peak in H2 S(0) intensity in the inner northwest spiral arm by 10". The H2 S(2) 
and H2 S(3) maps show H2 emission in the nucleus, spiral arms, and inter-arm regions of 
M51 and bar structure aligned north-to-south is apparent in H2 S(3) emission. H2 S(4) and 
H2 S(5) emission is resolved in the nuclear region of M51. 

2. The different distributions of H2 emission in M51 indicate significant spatial variations in 
H2 temperature and surface density. Using the low J (S(0) - S(2)) lines to trace the warm (T 
= 100 — 300 K) H2, we find that the warm H2 temperature is highest in the nuclear region 
at 192 K and the warm H2 surface density peaks in the northwest inner spiral arm at 11 Mq 
pc-2. Using the higher J (S(2) - S(5)) lines to trace the hot (T = 400 - 1000 K) H2, we find 
that the hot H2 temperature is lowest in the inner spiral arms (500 — 550 K) and increases to 
~ 600 K in the nucleus where the largest hot H2 surface density is found to be 0.24 Mq pc~^. 

3. The warm and the hot H2 surface density distributions are not cospatial and the warm- 
to- hot surface density ratio varies across M51. The warm H2 surface density distribution 
peaks in the northwest spiral arm and is offset from the hot mass peak by 11". The hot 
H2 surface density distribution shows two peaks, one in the nucleus of M51 and one located 
interior to the northwest inner spiral arm of M51. Variations in the warm-to- hot H2 ratio 
and differences in the distributions of the H2 line emission across M51 suggest that the warm 
H2 is mostly produced by UV photons in star forming regions while the hot H2 is mostly 
produced by shocks or X-rays associated with nuclear activity. 

4. The warm H2 follows the cold molecular gas traced by CO in the spiral arms of M51 
indicating that the warm phase is associated with the surface layers of dense molecular 
clouds. The H2 S(0) — H2 S(3) contours trace the CO; however, within the spiral arms, the 
peaks in H2 can be offset from the peaks in CO intensity. 

5. Comparing the distributions of H2 to V band imagery and H« emission reveals that 
the warm and hot H2 in the spiral arms is found in the dust lanes rather than spatially 
coincident with the Ha emission. 

6. [O IV] (25.89 /xm) emission and X-ray intensity peak in the nuclear region of M51 and 
their peaks are spatially coincident with the peak in hot H2 surface density. This implies 
that the hot H2 is more efficiently heated by the AGN, shocks (possibly associated with the 
AGN), or X-rays associated with the AGN. The spatial distributions of the [O IV] emission 
and X-ray surface brightness are very similar preventing the characterization of the primary 
excitation mechanism (shocks or X-rays) of the hot H2 phase in the nuclear region. 
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Fig. 1. — Overlay of the IRS footprints on the Spitzer IRAC 8 /xm image of M51. The three 
parallel strips mark the SL (5 - 14.5 /xm) observations with the white strip denoting the 
SL2 (5 - 7.5 iim) coverage and the grey strip denoting the SLl (7.5 - 14.5 //m) coverage. 
The central strip is where the SLl and SL2 overlap along with the LL (14 - 38 //m) radial 
strip giving us complete coverage of the mid-infrared spectrum across the nucleus and disk 
of M51. The small off-galaxy strip is the SL background observation. 
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Fig. 2. — Shown are sample Spitzer IRS low resolution spectra and the PAHFIT spectral 
decomposition for 3 different single pixel regions across the M51 strip. The SL and LL 
spectra are plotted separately because we decomposed the SL and LL data cubes separately. 
The locations of the regions are marked on the maps of the H2 surface density in Figure O 
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Fig. 3.— Maps of the H2 S(0) {top left), H2 S(l) {top right), H2 S(2) {middle left), H2 
S(3) {middle right), H2 S(4) {bottom left), and H2 S(5) {bottom right) line fluxes across 
the SL and LL strips that we mapped with the Spitzer IRS. The H2 S(0) and H2 S(l) maps 
are created from the LL data cube. The H2 S(2), H2 S(3), H2 S(4), and H2 S(5) maps are 
created from the SL data cube. The grayscale is in units of W m~^. Contour levels are at 
3.7 X 10"^^ 7.3 X 10-1^ 1.1 X 10-l^ 1.5 x 10-^^ 1.8 x 10'^^, 2.2 x 10-^^ and 2.9 x 
10-18 W m-2 for H2 S(0); 1.1 x IQ-^^ 2.1 x IQ-^^, 3.2 x IQ-^^, 4.3 x IQ-^^, 5.4 x IQ-^^ 
6.4 X IQ-i^ 7.5 X IQ-i^ 8.6 x IQ-^^ and 9.6 x IQ-^^ W m^^ ^2 S(l); 2.2 x IQ-^^, 4.4 
X 10-l^ 6.7 X 10-1^ 8.9 X 10-1^ and 1.1 x IQ-^s w m"^ for H2 S(2); 1.3 x IQ-^^, 4.0 x 
10-l^ 6.7 X 10-l^ 9.4 X 10-l^ and 1.2 x lO'^^ W m'^ ^2 S(3); 1.0 x lO'^^ and 2.0 
X 10-18 W m-2 for H2 S(4); 8.0 x 10"!^ 4.0 x lO'^^, and 7.3 x lO'^^W m'^ for H2 S(5). 
The vertical axis is the right ascension and the horizontal axis is the declination. In all of 
the maps, north is up, east is to the left, and the cross denotes the nucleus of the galaxy. 
The different spiral arm regions are labeled on the H2 S(0) and H2 S(l) maps in order to 
aid in discussing the molecular gas morphologies. The bar in the bottom right corner of the 
maps represents 1 kiloparsec. The box around the maps represents the SL or LL strip that 
we mapped. 
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Fig. 4. — Shown are excitation diagrams and the fits to the warm and hot H2 phases taken 
from three different single pixel regions along the M51 strip. The different regions are marked 
on the H2 surface density maps in Figure [51 
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Fig. 5.— Shown are the warm (T = 100 - 300 K) H2 {left) and hot (T = 400 - 1000 K) H2 
[right) surface density distributions. The surface density distributions are in units of Mq 
pc^^. Contours are overplotted for clarity. The warm H2 surface density contour levels are 
at 1.10, 2.21, 3.32, 4.43, 5.55, and 8.85 pc~^. The hot H2 contour levels are at intervals 
of 10 % of 0.25 M© pc^^. The hot H2 surface density distribution is derived from the fit to 
the H2 S(2) — H2 S(5) lines and the warm H2 surface density distribution is derived from 
the fit to the H2 S(0) — H2 S(2) lines, corrected for the contribution of the hot H2 phase. 
The three circles denote the regions of the spectra and excitation diagrams in Figures [2] and 
m respectively. 
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Fig. 6. — The warm (T = 100 — 300 K) H2 surface density distribution (in contours) 
compared to the warm H2 temperature distribution (in grayscale, in units Kelvin). The warm 
H2 temperature and surface density distributions are derived from the fit to the excitation 
diagrams across the strip for the H2 S(0) — H2 S(2) lines, corrected for the contribution of 
the hot (T = 400 — 1000 K) H2 phase. Surface density contour levels are at 1.10, 2.21, 
3.32, 4.43, 5.55, and 8.85 Mq pc~^ (same as in Figure E]). The non-rectangular shape of 
the map is due to the slight offset of the Spitzer IRS SL strip relative to the LL strip. In 
both the northwest and southeast corners of the maps, the temperature and surface density 
distributions could not be determined due to no detection of the H2 S(0) line. 
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Fig. 7. — The hot (T = 400 — 1000 K) H2 surface density distribution (in contours) com- 
pared to the hot H2 temperature distribution (in grayscale, in units Kelvin). The hot H2 
temperature and surface density distributions are derived from the fit to the excitation dia- 
grams across the strip for the H2 S(2) — H2 S(5) lines. Surface density contour levels are at 
intervals of 10% of 0.25 Mq pc~^ (same as in Figure ED- The high temperature and surface 
density of the hot phase towards the northwest and southeast ends of the strips is due to the 
lack of detection of extended emission from the H2 S(4) and H2 S(5) lines. In these regions 
(beyond the northwest and southeast spiral arms) we are unable to accurately determine the 
hot phase temperature and surface density. 
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Fig. 8. — The warm (T = 100 — 300 K) H2 surface density (in grayscale) compared to the 
hot (T = 400 — 1000 K) H2 surface density (in contours). Contours levels for the hot H2 
surface density distribution are at intervals of 10% of the maximum surface density (0.25 
M0 pc~^). The grayscale is in units of pc~^. 
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Fig. 9. — Left: Comparison of CO intensity (in grayscale) to the warm (T = 100 — 300 
K) H2 surface density (in contours). The CO intensity is in units of Jy km s~^. The warm 
H2 surface density contours are the same as in Figures [5] and El Right: Comparison of CO 
intensity (in grayscale) to the hot (T = 400 — 1000 K) H2 surface density (in contours). The 
CO intensity is in units of Jy km s~^. The hot H2 surface density contours are the same as 
in Figures [5] and [71 
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Fig. 10. — Comparison of the CO emission to the H2 S(0) {top left), H2 S(l) {top right), 
H2 S(2) {bottom left), and H2 S(3) {bottom right) emission. The CO emission maps are in 
units of Jy km s^^. Contour levels for H2 S(0), H2 S(l), H2 S(2), and H2 S(3) are the same 
as in Figure [31 



- 28 - 



15 



14 



13 



12 



11 



1 I 1 M • M I I t • •' ' I ' I ' ' I ' '" I" 



10 - 



47 09 





r 



I - 1 . i > 1 1 . 1 



13'^0'"05' 30"00' 55' 50^ 45' 40' 



0.3 



0,1 



15 



14 



13 



12 



11 



10 - 



47 09 



I ' ' ' ' I ' ' 




.1 




A. 



I . I 1 



0.2 



13''30°'05' 30™00° 55" SO' 45* 40' 



Fig. 11. — Left: Comparison of Ha (in grayscale) to the warm (T = 100 — 300 K) H2 
surface density (in contours). The Ha image is in units of counts s~^. The warm H2 surface 
density contours are the same as in Figures [5] and [61 Right: Comparison of Ha (in grayscale) 
to the hot (T = 400 — 1000 K) H2 surface density (in contours). The Ha image is in units 
of counts s~^ The hot H2 surface density contours are the same as in Figures [5] and [71 
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Fig. 12. — Left: Comparison of the V band image (in grayscale) to the warm (T = 100 — 
300 K) H2 surface density (in contours). The V band shows the dust lanes within M51. The 
V band image is in units of counts s~^. The warm H2 surface density contours are the same 
as in Figures and O Right: Comparison of V band image (in grayscale) to the hot (T = 
400 — 1000 K) H2 surface density (in contours). The V band image is in units of counts s~^. 
The hot H2 surface density contours are the same as in Figures [5] and [71 
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Fig. 13. — Comparison of Ha emission to the H2 S(0) {top left), H2 S(l) (top right), H2 
S(2) {bottom left), and H2 S(3) {bottom right) emission. The Ha image is in units of counts 
s"^ Contour levels for H2 S(0), H2 S(l), H2 S(2), and H2 S(3) are the same as in Figure [3l 
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Fig. 14. — Comparison of the V band image to the H2 S(0) {top left), H2 S(l) {top right), 
H2 S(2) {bottom left), and H2 S(3) {bottom right) emission. The V band image is in units 
of counts s~^. Contour levels for H2 S(0), H2 S(l), H2 S(2), and H2 S(3) are the same as in 
Figure [31 
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Fig. 15. — Left: Comparison of the [O IV] (25.89 fim) emission (in grayscale) to the warm 
(T = 100 K - 300 K) H2 surface density distribution (in contours). Hot H2 surface density 
contours are the same as in Figures and O Right: Comparison of the [O IV] (25.89 /im) 
emission (in grayscale) to the hot (T = 400 - 1000 K) H2 surface density distribution (in 
contours) . Hot H2 surface density contours are the same as in Figures [5] and [71 The [O 
IV] (25.89 fim) emission is in units of W m~^. 
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Fig. 16. — Left: Comparison of the smoothed 0.5 — 10 keV X-ray emission band (in 
grayscale) to the warm (T = 100 — 300 K) H2 surface density distribution (in contours). 
The X-ray image has been smoothed to the same resolution as the warm H2 surface density 
map. X-ray emission is in units of counts. The warm H2 surface density contours are the 
same as in Figures [5] and El Right: Comparison of the smoothed 0.5 — 10 keV X-ray 
emission band (in grayscale) to the hot (T = 400 — 1000 K) H2 surface density distribution 
(in contours). The hot H2 surface density distribution contours are the same as in Figures [5] 
and [3 
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Fig. 17. — Comparison of the 0.5 — 10 keV X-ray emission band (in grayscale) to the H2 S(2) 
{top left), H2 S(3) {top right), H2 S(4) {bottom left), and H2 S(5) {bottom right) emission 
in the nuclear region of M51. X-ray emission is in units of counts. The H2 S(2) and H2 S(3) 
emission contours are at 10 % of their peak values (2.20 x 10~^^ and 1.35 x 10~^^ W m~^, 
respectively). The H2 S(4) contours are at 1.0 x 10~^^ and 2.0 x 10~^® W m~^ and the H2 
S(5) contours are at 8.0 x 10"^^ 4.0 x IQ-^^, and 7.3 x lO'^^ W m-^. 



Table 1. H2 Parameters 



Transition 


Wavelength (/um) 


Rotational State (J) 


Energy (E/k) 


A (s-i) 


Statistical Weight (g) 


H2(0-0)S(0) 


28.22 


2 


510 


2.94x10"" 


5 


H2(0-0)S(1) 


17.04 


3 


1015 


4.76x10-1° 


21 


H2(0-0)S{2) 


12.28 


4 


1682 


2.76x10"^ 


9 


H2(0-0)S(3) 


9.66 


5 


2504 


9.84x10"'' 


33 


H2(0-0)S(4) 


8.03 


6 


3474 


2.64x10-8 


13 


H2(0-0)S(5) 


6.91 


7 


4586 


5.88x10-8 


45 



Note. — The statistical weight (g) is (2J +1)(2/+1) where / equals 1 for odd J transitions (ortho transitions) and I 
equals for even J transitions (para transitions). 



